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Abstract 

We calculate the in-medium D and Z)-meson masses in isospin aysmmetric nuclear matter in an 
effective chiral model. The D and D - mass modifications arising from their interactions with the 
nucleons and the scalar mesons in the effective hadronic model are seen to be appreciable at high 
densities and have a strong isospin dependence. These mass modifications can open the channels 
of the decay of the charmonium states Xa J/^) to DD pairs in the dense hadronic matter. 
The isospin asymmetry in the doublet D = (D°,D + ) is seen to be particularly appreciable at 
high densities and should show in observables like their production and flow in asymmetric heavy 
ion collisions in the compressed baryonic matter experiments in the future facility of FAIR, GSI. 
The results of the present work are compared to calculations of the D(D) in- medium masses in 
the literature using the QCD sum rule approach, quark meson coupling model, coupled channel 
approach as well as from the studies of quarkonium dissociation using heavy quark potentials from 
lattice QCD at finite temperatures. 
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I. INTRODUCTION 



The in-medium properties of hadrons is a topic of intense research in high energy physics 
since the last few decades. The ongoing and future relativistic heavy ion collision experi- 
ments, at the high energy accelerators SPS, CERN, Switzerland; SIS, GSI, Germany; RHIC, 
BNL, USA; LHC, CERN, Switzerland, and the compressed baryonic matter (CBM) ex- 
periments at the future facilities at GSI, Germany, are intended to probe matter at high 
temperatures and densities. The hadrons modified in the hot and dense hadronic medium 
resulting from heavy ion collision experiments, affect the experimental observables. E.g., the 
dilepton spectra observed from heavy ion collision experiments at the SPS are at- 

tributed to the medium modifications of vector mesons 

mum 

, and can not be 

explained by vacuum hadronic properties. The production of kaons and antikaons as well 
as their collective flow are attributed to the modifications of their spectral functions in the 

q, y, y, 0, fl q, 



medium 
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14j. Due to the presence of the light quark/antiquark, 



the medium modifications of the D (D) mesons can be appreciable in the dense hadronic 



medium 
collisions 



i^y h,q. 
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211 ] . The observation of open charm enhancement in nuclear 



22|,[23j as well as J/\l/ suppression as observed at the SPS 



24 



25 



261 ] . can be due 



to medium modifications of the D (D) in the medium. In higher energy heavy ion collision 
experiments at RHIC as well as at LHC, the J/\l/ suppression can arise from formation of 
a quark-gluon plasma (QGP) 27|, |28|]. However, the effect of the hadron absorption of J/^ 



is not negligible 29j, |30j, l3l|. With the drop in the mass of DD pair in the medium, the 



excited states of J/\l/, a major source of yield of J/\l/ 32], can decay to the D(D) final 



states 



33j, thereby leading to J/\l/ suppression 34| in the hadronic medium. The medium 



modification of the masses of J/\l/ and their excited states in the hadronic medium, due 



351 ] and the 



to the D-meson mass modifications have also been considered in the literature 
excited states of J^> are seen to have appreciable mass dop in the medium. The effects 
due to the mass modifications of the D-mesons as well as the charmonium bound states on 



experimental observables could be explored at the future accelerator facility at GSI 36]. 
It is thus important to study the modifications of the charmed mesons in the medium and 
hence one has to understand the charmed meson interactions in the hadronic phase. 

In the QCD sum rule approach, due to the presence of the light quark/antiquark, the 
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15,MQJ. In the 



2111. which is very similar to the value 



191 ]. Furthermore, lattice calculations 



mass modification of the .D-mesons arises from the light quark condensate 
quark meson coupling (QMC) model, the contribution from the m c (qq)N term is represented 
by a quark-cr meson coupling. The QMC model predicts the mass shift of the .D-meson to 
be of the order of 60 MeV at nuclear matter density 
obtained in the QCD sum rule calculations of Ref. 15 
for heavy quark potentials at finite temperature suggest a similar drop 

In this work we study the medium modification of the masses of open charm mesons 
(D 1 * 1 ) due to their interactions with the isoscalar-scalar mesons (a and ( mesons), isovector- 
scalar 5 mesons and the nucleons in asymmetric nuclear matter. The medium modification 
of the properties of kaons and antikaons in (isospin as ym metric) dense nuclear matter has 
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37]. 



38, 



391 ] and have been extended to 



been studied in a chiral S{7(3)-fiavor model in Ref. 
include the effects from hyperons in the asymmetric strange hadronic matter in ref. [40]. 
We generalise the model to S'[/(4)-flavor to derive the interactions of the D mesons with 
the nucleons and scalar mesons and investigate their mass modifications in the asymmetric 
matter. The masses of the D mesons in symmetric nuclear matter have been studied 



earlier in such an effective chiral model 



4l| . In a coupled channel approach for the study of 



D mesons, using a separable potential, it was shown that the resonance A c (2593) is generated 



dynamically in the 1=0 channel 



421 ] analogous to A (1405) in the coupled channel approach 



for the KN interaction 43J]. The approach has been generalized to study the spectral 
density of the D-mesons at finite temperatures and densitites 441 ]. taking into account the 
modifications of the nucleons in the medium. The results of this investigation seem to 
indicate a dominant increase in the width of the D-meson whereas there is only a very 



small change in the D-meson mass in the medium [44j . However, these calculations 42j, |44j , 
assume the interaction to be SU(3) symmetric in u,d,c quarks and ignore channels with 
charmed hadrons with strangeness. A coupled channel approach for the study of D-mesons 
has been developed based on SU(4) symmetry 45] to construct the effective interaction 
between pseudoscalar mesons in a 16-plet with baryons in 20-plet representation through 
exchange of vector mesons and with KSFR condition 0]. This model 45] has been modified 
in aspects like regularization method and has been used to study DN interactions in Ref. 



471 ] . This reproduces the resonance A c (2593) in the 1=0 channel and in addition generates 
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another resonance in the 1=1 channel at around 2770 MeV. These calculations have been 
generalized to finite temperatures accounting for the in-medium modifications of the 
nucleons in a Walecka type a — u model, to study the D and D properties [49] in the hot 
and dense hadronic matter. 

Within the effective chiral model considered in the present investigation, the D(D) en- 
ergies are modified due to a vectorial Weinberg- Tomozawa, due to scalar exchange terms 
(a, (,5) as well as range terms {39! I40]. The isospin asymmetric effects among D° and D + 
in the doublet, D= (D°,D + ) as well as between D° and D~ in the doublet, D = (D°,D~) 
arise due to the scalar-isovector 5 meson, due to asymmetric contributons in the Weinberg- 
Tomozawa term, as well as in the range term [3^]. 

We organize the paper as follows: We briefly recapitulate the S'L r (3_)-flavor chiral model 
adopted for the description of the asymmetric hadronic matter 39_|, |40j in Section II. The 
properties then are studied within this approach. These give rise to medium modifications for 
the D-masses through their interactions with the nucleons and scalar mesons as presented in 
Section III. Section IV discusses the results of the present investigation, while we summarise 
our findings and discuss possible outlook in Section V. 



II. THE HADRONIC CHIRAL SU(3) x SU(3) MODEL 



The effective hadronic chiral Lagrangian used in the present work is given as 

BW + £vec + A) + £-SB 

W=X,Y,V,A,u 



(1) 



are discussed. Eq. ([T]) corresponds to a relativistic quantum field theoretical model of 
baryons and mesons adopting a nonlinear rea 



broken scale invariance (for details see 



53 



54 



ization of chiral symmetry 



51 



52| and 



551 ]) as a description of the hadronic mat- 



ter. The model was used successfully to describe nuclear matter, finite nuclei, hypernuclei 
and neutron stars. The Lagrangian contains the baryon octet, the spin-0 and spin-1 meson 
multiplets as the elementary degrees of freedom. In Eq. (CD), Ckin is the kinetic energy 
term, Cbw contains the baryon-meson interactions in which the baryon-spin-0 meson inter- 
action terms generate the baryon masses. C vec describes the dynamical mass generation of 
the vector mesons via couplings to the scalar fields and contains additionally quartic self- 



interactions of the vector fields. Cq contains the meson-meson interaction terms inducing the 
spontaneous breaking of chiral symmetry as well as a scale invariance breaking logarithmic 
potential. Csb describes the explicit chiral symmetry breaking. 

The baryon-scalar meson interactions generate the baryon masses and the parameters 
corresponding to these interactions are adjusted so as to obtain the baryon masses as their 
experimentally measured vacuum values. For the baryon-vector meson interaction terms, 



there exist the F-type (antisymme 



the antisymmetric coupling 39, 



40 



; ric) and D-type (symmetric) couplings. Here we will use 
53 ] because, following the universality principle anc 



the vector meson dominance model, one can conclude that the symmetric coupling should be 
small. Additionally we choose the parameters 38, 
field d>, 



53J so as to decouple the strange vector 
sj^s from the nucleon, corresponding to an ideal mixing between u and (p. A 
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58 



59] has not been taken 



small deviation of the mixing angle from the ideal mixing 
into account in the present investigation. 

The Lagrangian densities corresponding to the interaction for the vector meson, C vec , the 
meson-meson interaction Co and that corresponding to the explicit chiral symmetry breaking 



38 



Csb have been described in detail in references 

To investigate the hadronic properties in the medium, we write the Lagrangian density 
within the chiral SU(3) model in the mean field approximation and determine the expectation 



values of the meson fields by minimizing the thermodynamical potential 54 



55]. 



III. D AND D MESONS IN THE MEDIUM 

We now examine the medium modifications for the D and D-meson masses in the asym- 
metric nuclear matter. The properties of nucleons and scalar mesons have been studied 
in the asymmetric hadronic matter within a chiral SU(3) model [39(. We assume that the 
additional effect of charmed particles in the medium leads to only marginal modifications 



60] of these hadronic properties and do not need to be taken into account here. However, to 
investigate the medium modification of the .D-meson mass, we need to know the interactions 
of the .D-mesons with the light hadron sector. 

The light quark condensate has been shown to play an important role for the shift in 
the D-meson mass in the QCD sum rule calculations [15]. In the present chiral model, 



the interactions to the scalar fields (nonstrange, a and strange, () as well as a vectorial 
Weinberg- Tomozawa interaction term modify the masses for D ± mesons in the medium. 
These interactions were considered within the SU(3) chiral model to investigate the modifi- 



39 



40]. 



cations of K-mesons in the dense (asymmetric) hadronic medium [38, 

To consider the medium effects on the D and D-meson masses we generalize the chiral 
S , f/(3)-flavor model to include the charmed mesons. The scalar meson multiplet has now 
the expectation value 



( (* + 5)/V2 






\ 



(2) 





(a - S)/y/2 

co 
v Q Q(cj 

with Q c corresponding to the cc-condensate. The pseudoscalar meson field P can be written, 
including the charmed mesons, as 



V2 



71 



7T — 

2K 



2K+ 2D° \ 
1+w 1+Wc 
2K° 2D~ 



y/2 1+w l+w c 
2K° 



1+w 1+w 
2D° 2D+ 













(3) 



/ 



\ 1 + Wc 1+Wc 

where w = \f2C,/<J and w c = ^/2( c /a. From PCAC, one gets the decay constants for the 
pseudoscalar mesons as f w = —a, fx = —{& + V2()/2 and fc = —(a + y/2(^ c )/2. In the 
present calculations, the value for the D-decay constant will be taken to be 135 MeV 19]. 
We note that for the decay constant of Df, the Particle Data Group 6j| quotes a value of 
f D + ~ 200MeV. Taking a similar value also for fjj would not affect our results qualitatively, 
however (see also js2]). 

The interaction Lagrangian modifying the D-meson mass can be written as 39J 



VN 



8fh 



+ (>7 M 



P 



3(jrfp + fi-fn) (D°(d,D°) - (d,D°)D°) + (D + (d,D-) - (d,D + )D- 
- n-fn) (D°(d,D ) - (d^D°)D°) - (D + {d,D-) - (8,D + )D- 



+ 



m 



D 



2f D 
1 



D 



[a + V2Q(D°D + {D-D+)) + 5(D°D°) - (D~ D+)) 



a + V2Q((d^D°)(d^D ) + (d„D-)(d»D- 



+ 8^(d ti D°)(d^D°) - (d ll D~)(d' J 'D + ) 



+ 



+ 



(pp + nn){{d ll D-)(d»D + ) + (d^D°)(d»D )) 
{pp + nn))({d^D°){d^D ) + (d f ,D-){d fi D + )) 



(pp — nn)((d tl D°)(d^D )) - (d ll D~)(d tl D + )) 



(4) 



In fll]) the first term is the vectorial interaction term obtained from the kinetic term in ([I]). 
The second term, which gives an attractive interaction for the D-mesons, is obtained from 
the explicit symmetry breaking term in (PQ) . The third term arises from kinetic term of the 

or the 



pesudoscalar mesons [39l . |4Cj. The fourth and fifth terms have been written down 



DN interactions, in analogy with the d\ and o?2 terms in chiral SU(3) model [39|, |40j. The 
last three terms in (j3J) represent the range term in the chiral model. It might be noted here 
that the interaction of the pseudoscalar mesons to the vector mesons, in addition to the 
pseudoscalar meson-nucleon vectorial interaction leads to a double counting in the linear 
realization of the chiral effective theory 63[]. Within the nonlinear realization of the chiral 
effective theories, such an interaction does not arise in the leading or sub-leading order, but 
only as a higher order contribution 631 ] . Hence the vector meson-pesudoscalar interaction 
will not be considered within the present investigation. 

The Fourier transformations of the equations of motion for D and D mesons yield the 
dispersion relations, 

to 2 + k 2 + m 2 D -U(uj, \k\) = 0, (5) 

where II denotes the self energy of the D (D) meson in the medium. 

Explicitly, the self energy Il(u;, \k\) for the D meson doublet, (D°,D + ) arising from the 
interaction (j3J) is given as 

1 



U(u,\k\) 



4/S 



Hpp + Pn) ± (pp- Pn)) 
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+ 



+ 



JD ^Jd 



d 2 



P% + P 3 n) ± (P S P - P 3 n, 



or 



k 



(6) 



7 



where the ± signs refer to the D° and D + respectively. In the above, &'{= o — ao), Cc'(= 
Cc — Ceo) an< i ^'( = $ ~ $o) are the fluctuations of the scalar-isoscalar fields a and £ c , and 
the third component of the scalar-isovector field, 5, from their vacuum expectation values. 
The vacuum expectation value of 5 is zero (<$o=0), since a nonzero value for it will break the 
isospin symmetry of the vacuum ( we neglect here the small isospin breaking effect arising 
from the mass and charge difference of the up and down quarks). In the above, p p and p n 
are the number densities of proton and neutron and p s p and p s n are their scalar densities. 
Similarly, for the D meson doublet, (D°,D~), the self-energy is calculated as 

1 



IL(u,\k\) 



'4/1, 



3 (Pp + Pn) ± (Pp ~ Pn) 



+ 
+ 
+ 



m 



D 



21 



D 



(a' + V2( c ' ± 5') 

JD ^J D 



((p S p + p S n)±(p 



S - P S 
r r 



(u 2 - k 



(7) 



where the ± signs refer to the D° and D respectively. 

The optical potentials are calculated from the energies of the D and D mesons 



U(u,k) = u(k) - \/k 2 + m 



2 



(8) 



where mo is the vacuum mass for the D(D) meson and u(k) is the momentum dependent 
energy of the D(D) meson. 



The parameter s di and c?2 are determined by a fit of 



scattering lengths [64l. 165. 



;he empirical values of the KN 



66| for 1=0 and 1=1 channels [39 



3- 



In the next section, we shall discuss the results for the D-meson mass modification ob- 
tained in the present effective chiral model as compared to the results in the literature, 
obtained from other approaches. 



IV. RESULTS AND DISCUSSIONS 

To study the D(D)-meson masses in asymmetric nuclear medium due to its interactions 
with the light hadrons, we have generalized the chiral SU(3) model used for the study 
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of the dense hadronic matter to SU(4) for the meson sector. The contributions from the 
various terms of the interaction Lagrangian (j4j) to the masses of the D = (D°, D + ) and 
D = (D°, D~) are shown in figures 1 and 2, as functions of density. These are illustrated for 
the isospin asymmetric case with the value of the asymmetry parameter, rj = (p n — p p ) / (2p#) 
as 0.5 and compared with the results obtained for the isospin symmetric case (r]=0). The 
isospin symmetric part (~ (p n + p p )) of the first term of (j3J), called the Weinberg- Tomozawa 
term, is attractive for D = (D°, D + ) mesons and leads to a drop of the masses of the D + 
and D° mesons, whereas it is repulsive for the D mesons in the nuclear medium, leads to 
an increase of the D° and D~ meson masses. In the isospin asymmetric nuclear medium, 
the Weinberg- Tomozawa term, leads to a mass splitting of the D° and D + mesons, giving 
a further drop in the mass of D + , whereas the asymmetry reduces the drop of the mass in 
D°. The D-meson self energy arising from the Weinberg- Tomozawa interaction, Hw T (u, \k\) 
is given by the first term of equation (jHJ), and at low densities, this turns out to be much 
smaller than (k 2 + m 2 D ). One can then, as a first approximation replace Hwt(&, \k\) by 
Uwri^D, \ k\) and solve for dispersion relation given by (jSJ). Confining our attention to the 
Weinberg- Tomozawa interaction only, the energies of the D° and D + mesons, in the above 
approximation, are given by 

u(\k\) ~ (\k\ 2 + ml) 1 ' 2 - -L [3(p p + Pn ) =F ( Pp - p n )} (9) 

One can note from the above equation (Q that at low densities, a given isospin asymmetry 
introduces equal increase (drop) for mass of D°(D + ) meson at low densities. However, at 
higher densities, there are deviations from the analytical expressions given by ([9]) as ex- 
pected, though the qualitative feature of the D + (D°) experiencing an attractive (repulsive) 
interaction from the vectorial Weinberg- Tomozawa term of (j6]) still remains the same, as 
can be seen from figure 1. One sees, from figure 2, that there is an increase in the masses 
of the D~ and D° and their mass shifts are equal for the case of isospin symmetric matter 
(r/=0). However, in the presence of asymmetry, the D~ mass is seen to have an increase 
whereas D° mass drops in the asymmetric medium. These behaviours for D~ and D° can 
be understood by examining the asymmetric contributions of the Weinberg- Tomozawa term 
(the first term of ((7j)). 

The scalar meson exchange contribution to the self energy is given by the second term of 

9 



equation ([6]) for D mesons and by the second term of (jTJ) for D mesons. Its interaction is 
attractive and is identical for both the D and D doublets for the isospin symmetric nuclear 
matter (a negligible difference in the energies is due to the difference in their vacuum masses, 
m,D+ = Wd-=1869 MeV and m^o = 771^0=1864.5 MeV). One might notice from the self 
energy terms due to the scalar meson exchange that a nonzero value of the scalar isovector 
5-meson arising due to isospin asymmetry in nuclear matter, gives a drop in the masses for 
the D + and D~ , whereas the 5 contribution is repulsive for D° and D° (a' = a — a > 
and 5' = 5 < 0). This gives the D° mass to be higher than D + mass for asymmetric nuclear 
matter as seen in figure 1 and the mass of D° mass (identical to D°) to be higher than D~ 
mass (identical to D + mass) as plotted in figure 2. However, one might notice that the shifts 
in the D + (D~) and D°(D°) masses about the isospin symmetric case {i]=0) case are not 
equal and opposite. The reason for the seen asymmetry in the mass splittings is the following. 
For low baryon densities, one has a' ~ p s ~ ps, and 5' ~ (p s p — p s n ) ~ (p p — p n ), so that 
one would expect the splittings of the masses of D + and D° to be symmetrical (about the 
isospin symmetric matter) for a given baryon density ps and isospin asymmetry parameter, 
77. However, these no longer hold good for higher densities and the mean field a calculated in 
the isospin symmetric situation (i]=0) turns out to be different for the asymmetric situation 
when the coupled equations of motions are solved for the scalar mean fields due to the 
presence of the scalar isovector 5 field, as compared to when these equations are solved for 
symmetric nuclear matter (for f]=0) in the absence of the 5 meson. 

The contributions to the D and D self energies due to range terms are given by the last 
three terms of the right hand side of equations ([6]) and ([7]) respectively. The first term of 
the range terms given by the third term in ([6]) is repulsive whereas the second and third 
range terms have attractive contributions, when the isospin asymmetry is not taken into 
account. However, the isospin asymmetry, due to a nonzero value of the 5 field, leads to an 
increase in the masses of the D + and D~ mesons and a drop in the masses of D° and D° 
from the isospin symmetric case. The second of the range terms (the d\ term) is attractive 
and gives identical mass drops for D + and D° in the D doublet as well as for D~ and D° 
in the D doublet. This term is proportional to (p^ + p"), which turns out to be different 
for the isospin asymmetric case as compared to the isospin symmetric nuclear matter, due 
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to the presence of the 5 meson. This is because the equations of motion for the scalar fields 
for the two situations (with/without 5 mesons) give different values for the mean field, a 
(~ (p s p + p s n ))- The last term of the range term (the d 2 term) has a negative contribution 
for the energies of D + and D° mesons as well as for D~ and D° mesons for the isospin 
symmetric matter. The isospin asymmetric part arising from the {{p s n — Ps P )) term of the 
di term has a further drop in the masses for D ± mesons, whereas it increases the masses 
of the D° and D° mesons from their isospin symmetric values. One sees from figures 1 and 
2 that for both D + and D~ mesons, the modification in the masses arising from the range 
terms, for the asymmetric matter (^=0.5), as compared to the isospin symmetric matter is 
negligible due to the increase from the first two range terms almost cancelling with the drop 
due to the d 2 term. For the D° and 5° mesons, the mass is seen to increase as compared to 
the isospin symmetric matter, due to the increase due to the second and third range terms 
dominating over the drop due to first range term. The values of the paramters d\ and d 2 



are fitted from the kaon-nucleon scattering lengths 38|, |39( to be 2.56/m^ and 0.73/rriK 
respectively [40| and it is seen that the d 2 term has a smaller contribution as compared to 
the di term. At high densities, these attractive di and d 2 terms dominate over the first range 
term (repulsive) and this leads to a decrease of the masses of the D and D mesons. There 
is seen to be a substantial drop of D-meson masses at high densities due to the inclusion of 
this range term. 

The density dependence of the masses of the D mesons at specific values of the isospin 
asymmetric parameter, 7], are shown in figure 3. The isospin asymmetry is seen to give a rise 
(drop) in the masses of the D° (D + ) as compared to their masses in the symmetric matter. 
For the isospin symmetric nuclear matter (rj=0), the drop in the mass of the D + at ps = Po 
is about 81 MeV from its vacuum value of 1869 MeV and D~ meson mass also has a drop of 
about 30 MeV, giving a mass splitting between the D + and D~ mesons as about 51 MeV. 
A similar drop of the D-meson mass is also predicted by the QMC model 2jJ as well as 
for the isospin averaged D(D) mass in the QCD sum rule approach in 151 ] . In the recent 



QCD sum rule calculations 



16] , the mass drop for D + D is about 50 MeV and the mass 



splitting is about 90 MeV at p = po, whereas for the isospin symmetric nuclear matter, we 
obtain these values as 110 MeV and 50 MeV respectively. The coupled channel calculations 
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44j show a dominant modification of the D-meson width and only a small change in the D 



meson mass in the medium. When the DN interaction is taken to be a Tomozawa Weinberg 
interaction supplemented by a scalar- isoscalar interaction |48|], the mass modification of the 
D mesons is seen to be around 10 MeV at po and about 50 MeV for a density of 2po. For 
the D~ meson, there is an increase in the mass by about 10 MeV and 30 MeV at densities 
of po an d 2p respectively 

In the present investigation, the isospin dependence is seen to be quite prominent for D° 
as compared to D + . The isospin asymmetry in the medium is seen to give an increase in 
the D° mass, whereas it gives a drop for the D + mass as compared to the r]=0 case. For 
the isospin asymmetry parameter, 77=0. 5, the D° mass is seen to rise by about 23 MeV and 
105 MeV for p# as po an d 5po respectively from their values of 1783 MeV and 1441 MeV in 
isospin symmetric matter. For the same value of asymmetry parameter, 7], the D + is seen 
to have a mass drop of about 18 MeV and 39 MeV for ps = po and 5po respectively, from its 
r)=0 values. This strong isospin dependence of the D mesons should show up in observables 
like their production as well as flow in asymmetric heavy ion collisions planned at the future 
facility at FAIR, GSI. Figure 4 shows the masses of the D~ and D°. It is seen that at high 
densities, both the D~ and D° are seen to have an increase in their masses in the asymmetric 
nuclear matter, as compared to the isospin symmetric case and this rise in the masses are 
seen to be similar for both D~ and D°. For example, at nuclear matter density, the masses 
of D~ and D° are 1838 MeV and 1840 MeV respectively for 77=0.5, which are very similar 
to the values of 1839 MeV and 1834.5 MeV for the isospin symmetric matter. For ps = 5po, 
the masses of D~ and D° are 1690 MeV and 1698 MeV respectively, which are higher by 
about 31 MeV and 43 MeV from the isospin symmetric case. One sees that the values of 
D~ and D° masses remain very similar at a given isospin asymmetry. But it is seen that the 
density effects on these masses are quite appreciable (a drop of about 30 MeV for nuclear 
matter density and of about 180 MeV for a density about five times nuclear matter density 
for 77=0.5) . 

Figures 5 and 6 show the optical potentials for the D and D doublets as functions of 
the momentum. These are shown for densities po an d 5po- The isospin dependence of the 
optical potentials is seen to be quite significant for high densities for the doublet (D + ,D°), 
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as has also been seen in the case of their masses. But, we see the optical potentials for 
the D have very similar values for D~ and D° for a fixed value of the isospin asymmetric 
parameter. But, as already has been seen for the case of their masses, we see the optical 
potentials for both D~ and D° are quite different from the symmetric nuclear matter case 
at high densities. 

The decay widths of the charmonium states can be modified by the level crossings between 
the excited states of J/\l/ (i.e., ^',Xc) and the threshold for DD creation due to the medium 
modifications of the D-meson masses [18[. In the vacuum, the resonances above the DD 
threshold, for example the state, has a width of 25 MeV due to the strong open charm 
channel. On the other hand, the resonances below the threshold have a narrow width of 
a few hundreds of KeV, only. With the medium modification of the D(D)-meson masses, 
the channels for the excited states of J/^f, like Xc ^' decaying to D + D~ or D°D° pairs 
can open up in the dense hadronic medium. This can increase the decay widths of \&' and 
Xc states at high densities. In figure 7, we show the density dependence of the mass of 
the D + D~ as well as D°D° pair, calculated in the present investigation. We see that the 
decay to D + D~ channel is almost insensitive to isospin asymmetry. On the other hand, the 
decay channel to the D°D° is seen to have a strong isospin asymmetry dependence, with the 
asymmetry shifting the onset of the decay to higher densities. The decay of charmonium 

331 ] . It is seen to depend sensitively on the relative 



states to DD has been studied in Ref. 



18 



momentum in the final state. These excited states might become narrow [18( though the 
D-meson mass is decreased appreciably at high densities. It may even vanish at certain 
momenta corresponding to nodes in the wavefunction [isj]. Though the decay widths for 
these excited states can be modified by their wave functions, the partial decay width of Xc2, 
due to absence of any nodes, can increase monotonically with the drop of the D + D~ pair 



mass in the medium 



181 ]. This can give rise to depletion in the J/ty yield in heavy ion 



collisions. The dissociation of the quarkonium states (^', Xc J /^) into DD pairs have also 
),0] by comparing their binding energies with lattice results on temperature 
dependence of heavy quark effective potential [3j 



been studied 120 
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V. SUMMARY 



To summarize we have investigated in a chiral model the in-medium masses of the D, D- 
mesons in asymmetric nuclear matter, arising due to their interactions with the nucleons 
and scalar mesons. The properties of the light hadrons - as studied in a SU(3) chiral model 
- modify the D(D)-meson properties in the dense hadronic medium. The SU(3) model with 
paramaters fixed from the properties of hadron masses in the vacuum, and low energy KN 
scattering data, is extended to SU(4) to derive the interactions of D(D)-mesons with the 
light hadron sector. The mass modifications for the D mesons are seen to be similar to 
earlier finite density calculations of QCD sum rules [if], Q] as well as to the quark meson 
coupling (QMC) model 21], in contrast to the small mass modifications in the coupled 



channel approach 



44 



48] . In our calculations, the presence of the repulsive range term (the 



fourth term of (j4]) is compensated by the attractive d\ and di terms given by the last two 
terms in (j4]) and the latter (dominated by the d\ term) have an effect of reducing the masses 
of both D and D mesons at high densities. 

The medium modification of the D-masses can lead to a suppression in the J/^f- yield in 
heavy- ion collisions, since the excited states of J/\l/ and at a much higher density (~ 5po), 
J/ty can decay to DD pairs in the dense hadronic medium. The decay to D + D~ pair seems 
to be insensitive to isospin dependence, whereas due to increase in the mass of D°D° in the 
asymmetric medium, isospin asymmetry is seen to disfavour the decay of the charmonium 
states to D°D° pair. The isospin dependence of D + and D° masses is seen to be dominant 
medium effect at high densities which might show in their production (D + /D°), whereas for 
the D~ and D°, one does see that even though these have a strong density dependence, their 
in-medium masses remain similar at a given value for the isospin asymmetry parameter t]. 
The strong density dependence as well as the isospin dependence of the D(Z)) meson optical 
potentials in the asymmetric nuclear matter can be tested in the asymmetric heavy ion 
collision experiments at the future GSI facility [36]. 
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